A solution in 3 dimensions for current in a semiconductor under high level injection from a point 

contact 



Mike W. Denhoff* 

National Research Council, Institute for Microstructural Sciences, 
M-50, Montreal Road, Ottawa, ON, K1A 0R6, Canada 
(Published in: 28th Intl. Conf. on the Physics of Semiconductors, Vienna, 24-28 July 2006, AIP Conference Proceedings Vol. 893.) 

The standard equations for semiconductor device analysis were solved by specifying the electron and hole 
current injected at a small contact, assuming high-level injection. Calculated current-voltage characteristics 
were fit to measurements of a single point breakdown in an ultrathin dielectric. It was found that the minority 
carrier injection level was about 70%. 
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INTRODUCTION 

The nature of the contact to a semiconductor at point break- 
downs in ultrathin oxides, is not well understood. An impor- 
tant observed feature of these contacts is that they can have a 
high minority carrier injection ratio. There have been a num- 
ber of attempts to model point breakdowns in ultrathin oxides 
on silicon [1]. These models concentrate on the physics at 
the contact interface and attempt to derive the current vs volt- 
age (IV) characteristic. However, in many physical devices 
the main contribution to the IV characteristic is due to the 
conduction in the bulk of the semiconductor. Therefore, one 
should include a model of the bulk conductivity in an expla- 
nation of point contact behavior. In this paper, such a model is 
developed using standard semiconductor equations solved for 
point contacts with minority carrier injection. 



POINT CONTACT PHYSICS 

This problem has the same geometry as the spreading resis- 
tance problem. The difference is that the spreading resistance 
problem assumes an ohmic contact and only majority carrier 
conduction. In the point contact problem arbitrary levels of 
minority carrier injection will be included. This calculation 
will assume high-level injection, which implies that there is 
charge neutrality in the bulk. This means that the number of 
excess holes (beyond the equilibrium concentration) is equal 
to the number of excess electrons. Since the bulk is neutral, 
the Laplace equation can be used to find the electric field, 



VE = 0. 



(1) 



The remaining equations needed to describe conduction in a 
semiconductor are the continuity equations for holes and elec- 
trons along with expressions for the hole and electron current 
in terms of drift and diffusion [2, 3]. For high-level injection, 
as in the case of high forward bias of a pn junction [2], these 
equations reduce to 
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for holes and electrons. p t (n t ) is the total hole (electron) con- 
centration, p (n) is the excess hole (electron) concentration, 
and po (no) is the equilibrium hole (electron) concentration, 
so that pr — p + po (n t = n+ no)- L p is the high-level diffu- 
sion length, Li = 2jj,„D p x p /(jJ,„ + /i p ), where fi p and \i n are 
hole and electron mobilities, D p is the hole diffusion constant, 
and % p is the hole lifetime. A similar equation can be written 
for L„. Since ji p D„ — ji n D p and x p = T„, then L p = L„ = L. 

These equations must be solved for the boundary conditions 
of the contact. A semi-infinite structure will be used, where 
the surface of the semiconductor is on the x-y plane and the 
bulk of the semiconductor extends in the positive z direction to 
infinity. The potential at infinity will be taken as 0. There are 
mixed boundary conditions on the x-y plane. Similar mixed 
boundary conditions, for the Laplace equation, are used in the 
solution of the usual spreading resistance problem [4, 5]. Out- 
side the contact radius a, the current through the surface and 
the perpendicular component of the electric field are 0. The 
boundary conditions on the contact are given by specifying the 
hole and electron currents through the contact. These condi- 
tions are set on the semiconductor side of the contact in order 
to avoid the physics of the contact interface. 



3-D SOLUTION 

The problem is cylindrically symmetric about the z-axis 
which is perpendicular to the center of the contact. Equa- 
tion (2) can be written (ignoring the 9 dependence) in cylin- 
drical coordinates, 
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This is the Helmholtz equation. It can be solved by separating 
variables. Using the boundary condition that the potential is 
zero for r, z at infinity, the solution for excess hole concentra- 
tion is 

p(r,z) = [ A(X)J (Xr)exp(-[X 2 +l/L 2 ]?z)dX. (4) 
Jo 

Jo is a Bessel function of order 0. The function A(X) must 
be found using the remaining boundary conditions. A similar 
solution for the electron density can be written. 
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single point breakdown data 
a=1 nm, 75% electrons 
a=10 nm, 74.33% electrons 
a=50 nm, 100% holes 
a=50 nm, 71% electrons 
a=50 nm, 71.5% electrons 
a=50 nm, 72% electrons 
a=50 nm, 80% electrons 
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FIG. 1: IV data of a single point breakdown. Theory curves have 
voltage offsets of -2.8 V, -4.5 V, and -5.2 V for a = 1 nm, a = 10 nm, 
and a = 50 nm, respectively. 



An analytic solution can be found if A (A) is chosen to be 



A{X) ~ B (A2 + C 2) A ' 



(5) 



where J\ is a Bessel function of order 1, B is a constant, and 
c 2 = 1/L 2 . With this choice, the current density over the area 
of the contact is approximately constant for aCI. The elec- 
tric field at the contact must be consistent with the current, 
which implies that the on the area of the contact the z compo- 
nent of the field is 



E z = EJo{c\fa?- 



-r 2 ), for r < a 
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where £ is a constant. Since, in the high-level injection case, 
p = n, the solutions can be equated and values of B and E 
determined for any given values of total hole current, //,, and 
electron current, /„. The potential at the center of the contact 
can be found from the solution of the Laplace equation and 
thus the / vs V curve can be calculated. The numerical results 
below use typical values of parameters for silicon. 



Single point breakdown 

A single point breakdown device was fabricated by open- 
ing a 100 nm square window in a thick (20 nm) SiC»2 insu- 
lating layer on a 6 £2-cm p-type silicon wafer. A thin SiC>2 
layer (2 nm) was grown in this opening and an indium tin ox- 
ide top electrode was deposited on top of this [6]. The device 
was voltage stressed to create a hard breakdown. IV data of 
this breakdown site is presented in Figure 1 . Theoretical IV 
curves were calculated assuming a constant electron injection 
ratio. A constant voltage was added to model the voltage drop 
across the contact interface for large current. The experimen- 
tal current is more than 10 times larger than the calculated 
current for only hole injection. Considering the curvature at 
low currents, the best fit is for a = 10 nm and an electron in- 
jection of 74%. With smaller values of a the theoretical curve 
is not straight and cannot fit the measured data. A reason- 
able fit can be obtained for a range of contact size, however 
it is clear that the electron injection ratio is between 70% and 
75%. This demonstrates that a large minority injection ratio 
exists for these devices. 
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Injection of electrons 

An interesting case is that of injecting a current of pure elec- 
trons into a p-type silicon. In the case of large injected current, 
the solution for E and the resistance can be approximated as 

„ -D„ 0.85D„ 

k = tt^ — r- y'> 
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G is somewhat dependent on a, but is approximately equal to 
4/(37r). It is seen that as the electron current gets large, the 
electric field coefficient tends to a constant. This is opposed 
to the normal Ohm's law behavior where the electric field is 
proportional to current. The resistance of the contact actually 
decreases as l/I„ as opposed to the ohmic behavior of a con- 
stant resistance. 



